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I. INTRODUCTION 



Two heavy antiquarks in a doubly heavy antibaryon feel an attractive force when they 
are in the 3 of color SU(3) and therefore are expected to form a compact diquark whose 
size, tqq ~ (mgf) -1 , is much smaller than Aq^ d . Here niq is the heavy quark mass and 
v is the typical velocity of the antiquarks in the diquark. To the light degrees of freedom 
in the doubly heavy antibaryon, this diquark is indistinguishable from a heavy quark in the 
rriQ — > oo limit. Therefore, the properties of ground state doubly heavy baryons should 
be related to properties of singly heavy mesons. Savage and Wise [l[ first observed that 
quark-diquark symmetry is a consequence of the heavy quark limit and used Heavy Quark 
Effective Theory (HQET) to derive a relation between hyperfine splittings of doubly heavy 
baryons and singly heavy mesons. 1 The correct framework for analyzing the properties of 
hadrons with two or more heavy quarks is Non-Relativistic QCD (NRQCD) [1, H, 0] and 
recently the heavy quark-diquark symmetry prediction was rederived in both the potential 
NRQCD Q and velocity NRQCD E l fo rmalisms. An extension of heavy hadron chiral 
perturbation theory (HHyPT) jloL I 111 . Il2j ] that incorporates heavy quark-diquark symmetry 
was developed in Ref. (l3( and used to derive a novel quark-diquark symmetry prediction 
for the electromagnetic decays of the spin-3/2 members of the ground state doubly heavy 
baryon doublets, as well as study 0(l/mg) and chiral corrections to the electromagnetic 
decays and doubly heavy baryon masses. Some quark-diquark symmetry predictions for 
the exotic heavy tetraquarks as well as a critical discussion of the applicability of heavy 
quark-diquark symmetry to charm hadrons appeared in Ref. Q. 



Much of the recent theoretical work has been stimulated by the SELEX collaboration 
reports of candidates for doubly charm baryons 0, r| 3| • While the masses and hyperfine 



splittings of the observed states are consistent with theoretical calculations in the quark 
model |4| as well as quenched lattice QCD, other aspects of the data are poorly understood at 
the present time. For example, isospin splittings, weak decay rates, and production yields are 
in disagreement with theoretical expectations. Some excited states which would be expected 
to decay electromagnetically are instead observed via their weak decay. For a review of the 
experimental situation and discussion of theoretical interpretations of the SELEX results, 
see Refs. O Q El] Recently, the BABAR Q and BELLE experiments searched for 



doubly charmed baryons in e + e~ collisions but failed to observe the SELEX states. More 
experimental data and improved theoretical understanding of doubly charm baryons are 
necessary to clarify the situation. 

Lattice gauge theory studies of the properties of doubly heavy baryons with dynamical 
quarks would be quite useful. First principles calculation of the spectroscopy of doubly 
heavy baryons could be directly compared with experiment to test our understanding of 
QCD and aid in the interpretation of candidate states. It would also be of interest to have 
first principle calculations of electromagnetic and weak decay matrix elements. Finally, 
the lattice could help determine the suitability of quark-diquark symmetry for charm and 
bottom hadrons. In these cases, it is possible that symmetry breaking corrections could be 
large enough to spoil symmetry predictions, but no systematic study of symmetry breaking 
corrections in NRQCD is yet available. A theoretical testing ground in which the heavy 
quark mass can be dialed at will could determine whether this symmetry will turn out to be 
useful for charm and bottom. Spectroscopic studies of doubly charm and bottom baryons 



1 The quark model has also been used to derive the relation 
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in quenched lattice QCD appeared in Refs. (3, Hi, 22|. The heavy quark-diquark symmetry 
prediction for the hyperfine splitting holds at the 20% level, which is reasonable given the 
sources of error. In QCD, corrections are expected to be 0{y 2 ) and/or 0(Aqcd/^q)- In 
existing lattice calculations there are also uncontrolled errors due to quenching, which, as we 
will see below, could be significant. Hopefully, the near future will bring studies of doubly 
heavy baryons in simulations with dynamical quarks. 

As a step towards such studies, it is clearly useful to develop extensions of the heavy quark- 
diquark symmetric version of HHxPT |l3| that include lattice artifacts such as quenching 
and partial quenching. These theories provide the formulae needed for chiral extrapolations 
in lattice calculations. We will focus on chiral corrections to the doubly heavy baryon masses, 
especially hyperfine splittings and the quark-diquark symmetry prediction. We also apply 
the formalism to the electromagnetic decays of the spin-3/2 members of the ground state 
doubly heavy baryon doublets. 

The organization of our paper is as follows. In section El we review the basics of 
quenched chiral perturbation theory (Q;\PT) and partially quenched chiral perturbation 
theory (PQ;\PT) in the pseudoscalar meson sector. We then write down versions of PQ%PT 
and QxPT that incorporate heavy quark-diquark symmetry. In section ITTTl we calculate the 
one-loop chiral corrections to the masses of doubly heavy baryons in PQxPT and QxPT. 
These expressions can be used to perform the chiral extrapolation of lattice data on the 
spectrum of doubly heavy baryons provided one is in the range of applicability of the chiral 
effective theory. In section ITY] we calculate to one loop the dipole and quadrupole transition 
moments in PQ%PT. For these transitions, the lattice can be used to test the predictions 
for the electromagnetic decay widths of the spin-3/2 members of the ground state doubly 
heavy baryon doublets ^3( . Finally, we conclude with a brief summary in section Some 
useful formuale are collected in the Appendix. 

II. PQxPT AND QxPT FOR DOUBLY HEAVY BARYONS 

The formalism for constructing effective theories for quenched and partially quenched 
QCD is now well established. We use the technique of graded Lie algebras that was for- 
mulated in Refs. @, 0, H H B HI Hj- Equivalently one could use the replica 



method [3l|. We begin by briefly reviewing the pseudoscalar meson sector of Q;\PT and 
PQxPT, and then write down the quenched and partially quenched chiral Lagrangians in- 
corporating heavy quark-diquark symmetry. These theories are an extension of QxPT and 
PQxPT for singly heavy mesons [32j - 

A. Pseudoscalar Mesons 

In the pseudoscalar meson sector, we can treat PQxPT and Q%PT as different limits of 
the same theory. This theory is described by the Lagrangian 

C = Il s tr (d»E%E) + ~ str (m q tf + mjs) + a*0"$ o fy$o - /zg$g, (1) 
o 4 



where the field E is defined by 

v (' 

x = 1 > V I 1 I - 



S = exp ( 2 -^\ = e, (2) 
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and the meson fields appear in the i7(6|3) matrix, 

The quantities a$ and /io are non- vanishing in the chiral limit. The M and M matrices 
contain bosonic mesons (with quantum numbers of qq pairs and qq pairs, respectively, where 
q is a valence or sea quark and q is a ghost quark), while the \ and x matrices contain 
fermionic mesons (with quantum numbers of qq pairs and qq pairs, respectively). The upper 
3x3 block of the matrix M contains the familiar pions, kaons, and eta, while the remaining 
components consist of mesons formed from one or two sea quarks. The operation str() in 
Eq. (JTJ) is a supertrace over flavor indices, i.e., str(A) = ^2 a e a A aa , where e a = (—) 1+Va , and 
r] a are the grading factors of the Lie algebra. 

The quark mass matrix appearing above is given by 

m q = diag(m„, m d , m s , uij, m h m r , m u , m d , m s ). (4) 

We will work in the isospin limit of both the valence and sea sectors, where = m u and 
mi = rrij. Expanding the Lagrangian in Eq. (JTJ) to lowest order in the fields, one finds that 
mesons with quark content qq' are canonically normalized and their masses are given by 

™< 2 qq > = jk( m q + m q')- ( 5 ) 

The propagators of the off-diagonal mesons have the usual Klein-Gordon form. 

On the flavor diagonal, the situation is more complicated and differs dramatically between 
the quenched and partially quenched theories. The flavor singlet field that appears above is 
defined to be $o = str(<&)/\/6. Using this, the leading-order r] a r]b propagator is 29} 

g , 2x = 1 ~ m )j) ~ m rr) ( a ^ 2 ~ /jj / g N 

ab Q 2 - rnl a + ie 3 (q 2 - m 2 aa + ie) (q 2 - m 2 bb + it) D(q 2 ) ' 

The function D(q 2 ) is defined as 

D(q 2 ) = {q 2 ~ m-,)(<? 2 - m 2 rr ) + (a%q 2 - [x 2 )(q 2 - m 2 x ), 

and the mass mx is given by m 2 x = | (jn 2 j + 2m 2 r ) . 

Because the fermionic determinant in quenched QCD is set to a constant, the sea quark 
masses are effectively infinite. The mesons containing sea quarks decouple from the theory, 
leaving a Z7(3|3) Goldstone manifold. In this limit, the flavor neutral propagator in Eq. © 
becomes 

X 2 2 
gQ fg2\ _ le a°ab — /i Q ^ 

ab q 2 - m 2 aa + ie 3 (q 2 - m\ a + ie) (q 2 -m 2 bb + ieY 

where we have appended a superscript Q to denote that this propagator is for the quenched 
theory. Notice the flavor singlet meson does not decouple as the theory still depends on the 
parameters a$ and /j, . Numerical data suggest that a$ is small (for example, a recent study 
finds ct$ = 0.03 ± 0.03 |33 | 1 , and we neglect it below; /i will be treated in the standard 
quenched power counting [34]. 
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In the partially quenched theory, sea quark contributions are retained and partially 
quenched QCD (PQQCD) has a strong axial anomaly. Because of this, the mass of the 
singlet field fio can be taken to be on the order of the chiral symmetry breaking scale, and 
the <3> integrated out of the theory. The Goldstone manifold is reduced to S77(6|3) but the 
resulting flavor neutral two-point function, however, still deviates from the familiar form in 
%PT. In PQQCD this propagator is given by the /iq — > oo limit of Eq. JBJ), namely 



a2V 



it J, 



ab 



(q 2 - m 2 ^) (q 2 - 



nr. 



<? 2 ~ m la + i€ 3 (? 2 ~ m la + i€ ) ((I 2 ~ m bb + ie ) ((I 2 -™ 2 X + le ) 

The PQxPT flavor neutral propagator can be conveniently rewritten as 

Gab = e a$abPa + Pab {Pa, Pbi Px) , 



(8) 



(9) 



where 



Tab (A,B,C) 



Q ~ m 2 aa + i e 



9 2 i • ' 

Q ~ m bb + ie 



(m 2 aa - m 2 j3 ) {m 2 aa - m. 



A + 



(mi 



rn 



[mi 



m 



xj 



{m\ - m%) 



(m 



x 



X 



9 2 i • ' 

q z — m x + ie 
(m 2 bb - m%) {m 2 bb - m 



[m 



bb 



mi 



[m 



bb 



m 



X) 



[m\ - m 2 aa ) 



C 



[m 



x 



m 



bb) 



B 



(10) 



B. Doubly heavy baryons 

As discussed earlier, in the heavy quark limit the properties of singlv heavy mesons and 
doubly heavy baryons are related by a U(5) quark-diquark symmetry Recently a chiral 
Lagrangian incorporating this U(5) symmetry has been derived flli ]. It is written in terms 
of a super-field TC given by 

7~La,nf3 = Ha,af3 + T a>i p , (11) 

where a is an anti-fundamental flavor index, a and (3 are Pauli spinor indices, while i is a 
vector spin index. The super-spin index fi = 1-5 and is related to a and i via 

t* = + + (i + ml + % + Sl). (12) 

The field H a ^ a/ 3 is the heavy meson field of HHxPT [HI El 2| i n the rest frame of the heavy 
meson: 

H a ,af3 = Pa ' °a/3 + PJ>a.fi- (13) 

H a contains the vector (P*) and pseudoscalar (P) field components. The field T a ^ is the 
doubly heavy baryon field, 

Ta,ip = V2 \ + -^Ha^o-^J , (14) 

that contains the spin-1/2 (H) and spin-3/2 (H*) field components. The spin-3/2 field 
satisfies a non-relativistic Rarita-Schwinger constraint, 2* i/3 o"^ 7 = 0. 
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In terms of the super-field 7i, the Lagrangian for the quenched (Q) theory is 
C Q = (H\WD Q ) \ - g Q (H j HA ■ a) - g' Q (H j H(r) ■ str(A) 



4 

while that for the partially quenched theory (PQ) is 
C p Q = («t(Wi5 )) -g p Q{rfHA-o) 



(^S • Ho-) + a Q (tfHM) + a' Q (V)U) str(A4), (15) 



a pq 

(H^ ■ Ha) + a PQ (H^HM) + a' PQ (H^H) MM). (16) 

Both Lagrangians include the leading quark-diquark symmetry breaking term. The bracket 
() notation denotes the contraction of flavor, spinor and super-spin indices of field bilinears. 
For a generic super-spin matrix X, flavor matrix Y, and spin matrix Z, () denotes 

{rCxuYz) = Y J Y,Y, e - n i^x^n b , u pY ba Zp a . (17) 

a,P a,b 

In Eqs. ()15p and (|16|). A is the spatial part of the axial- vector pion field, = 
2 ~ an< ^ M is the mass operator defined by M. = | (£,m q (, + ^Tn q ^). The 

U(5i) super-spin matrix I] is given by 



where (T 1 )^ = -ie^. 

Notice that there are separate Lagrangians for quenched and partially quenched QCD. 
The low-energy constants of the two theories are different, e.g. g p Q ^ gQ. Notice also that 
in QxPT the flavor singlet field cannot be integrated out. Consequently str(A) 7^ and the 
axial hairpin coupling g'Q remains. Additionally in Q%PT, str(A^) = to leading order, 
and so we can effectively set a'Q = for our calculations. There is no sector of QxPT that 
contains j^PT, so in general quenched observables are unrelated to QCD observables. The 
sea sector of PQQCD is QCD, and hence PQxPT contains the low-energy constants of x?T. 
Restricting the flavor indices to the sea sector in Eq. 1)160. we recover the chiral Lagrangian 
of Ref. [r|, so g p Q, A^ Q , o p Q and o' p Q are identical to the low-energy constants g, A H , a, 
and a' of xPT, respectively. At this order there are no additional operators in the partially 
quenched Lagrangian. Carrying out the super-spin sums in Eqs. (1150 and (|16|). we produce 
separate terms for heavy mesons and doubly heavy baryons but with the same low-energy 
constants. The heavy meson sector of the PQ^PT and QxPT Lagrangians above is that of 
Ref. 1321. 



III. MASSES 



In this section, we use the quenched and partially quenched chiral Lagrangians that 
incorporate quark-diquark symmetry to calculate chiral corrections to the masses of doubly 
heavy baryons. The results can be used for the chiral extrapolation of lattice data. PQ^PT 
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is considered in subsection III. A and Qx?T is discussed in subsection III.B. It is important 
to keep in mind that while PQxPT actually contains xPT, there is no such connection 
between Q%PT and x?T. 

Near the chiral limit, the masses of doubly heavy baryons generically have the form 

M s = M - ^A H + il4 1} + Mi 3/2) + . . . 

M H * = M + ^A H + M$ + M§J 2) + . . . , 

where M is the mass of the doubly heavy baryons in the chiral and heavy quark limit and 
Ajj is the hyperfme mass splitting of the heavy mesons in the chiral limit. 2 The hyperfme 
splitting of the H* and 5 baryons is |A# in the chiral limit due to quark-diquark symmetry. 
The term M^ 1 ' is proportional to m q and comes from tree-level counterterms. The term 
comes from the one-loop diagrams and depends on m, the Goldstone boson mass, 
and S, the mass splitting between the external and virtual doubly heavy baryons. Goldstone 
boson masses are oc y/m q , while 5 receives contributions from SU(3) breaking counterterms, 
which are oc m q , and/or hyperfme splittings, which scale as I/tuq. In our power counting 
we take 5 ~ m ~ Q in which case M^ 3 / 2 ) oc Q 3 . In the chiral and heavy quark limits, S can 
be neglected relative to m and then M^ 3 / 2 ) oc vrC'J 2 . The ellipsis denotes higher-order terms 
in the Q expansion. 



A. PQxPT 

At linear order in the quark mass, the contributions to the masses of doubly heavy baryons 
at tree level come from the mass operators in Eq. 

= Af£ = -a(m q ) a - a'stv(m q ), (19) 

where str(m 9 ) = 2rrij + m r in the isospin limit, and the subscript a denotes the light-quark 
anti-fundamental index of the doubly heavy baryon. Combining these results with the mass 
splitting in the chiral limit, we then define the baryon mass splittings 

A ba = M Sb - M Ea = M a » - M ss = a[(m q ) a - (m q ) b \, (20) 
A* ba = M sl - M Ho = ^A H + a[{m q ) a - {m q ) b \. (21) 

At one-loop order, the doubly heavy baryons receive mass contributions in PQxPT from 
the sunset diagrams shown in Fig. The vertices in these diagrams are generated from 



2 Strictly speaking Ah is not same the parameter appearing in the Lagrangian above, Eq. (|16J) . which in 
this footnote we denote by A^. This parameter, which transforms as a chiral singlet, is corrected by a 
string of operators that differ only by scaling with powers of A° H /A X . Thus we subsume this dependence 
as Ajj = Ajj(A° h ) including also contributions to Ajj arising from loop graphs. In this way we obtain 
Ajj as the chiral limit value. The same is true of all other low-energy constants in Eq. H16fl . We similarly 
drop their polynomial dependence on A^/A x , and absorb contributions from loop graphs to work with 
the values of the low-energy constants in the chiral limit. For more details see, e.g. Ref j^. 
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FIG. 1: One- loop graphs which give contributions to the masses of the doubly heavy baryons in 
PQxPT. The single and double lines correspond to 3 baryons or H* baryons, respectively, while 
the dashed lines correspond to mesons. The filled squares denote the non-singlet axial coupling 
given in Eq. Q16JI . while crosses denote the hairpin interaction. 

the pion-baryon interaction terms in Eq. (|16jh while the hairpin interaction corresponds to 
the disconnected part of the flavor neutral propagator in Eq. (JSJ). Calculating these PQ^PT 
loop diagrams, we find 



Mr 



(3/2) 



(4tt/)s 



b=j,l,r 



5 4 

-K(A ba , m ba , n) + -K(-A* ab , m ba , //) 



(4tt/) 2 

for the spin-3/2 doubly heavy baryons, H*, and 



5 4 
-/C(0, m aa , m aa , fj) + -JC[-A* aa , m aa , m aa , //) 



M (3/2) 



(4tt/) 5 



E 



b=j,l,r 



-K(A ba , m ba , fi) + -K{A* ba , m ba , p) 



(4tt/) 5 



-/C(0, m aa , m aa , //) + ~/C(A* a , m aa , m aa , //) 



for the spin-1/2 doubly heavy baryons, S Q . The non-analytic function K(S,m,n) is 



mi? 



5 



K{8,m,n) = 2(5 2 -m 2 ) 
where R(x) is defined by 

= V x 2 — 1 log 



5 log 



2\ 1 



+ 5m 2 log 



x — 



1+ie 



,i + \lx 2 — 1 + ie 

while the function K.(S,m aa ,m bb , /x) which arises from hairpins is given by 

/C(5, m aa , m 66 , /i) = P ab m aa , //), m 66 , //), m x , //) 
with P afe appearing in Eq. (jTDjl. and a limit is understood when b = a. 



(22) 



(23) 



(24) 



(25) 



(26) 



S 




m [ MeV ] 

7U sea 

FIG. 2: 5\hf as a function of m nsea , for different values of m-KvaX- The width of the bands is the 
result of varying fi between 500 MeV and 1500 MeV. 

To study the chiral non-analytic behavior of doubly heavy baryon mass splittings in 
PQxPT, we need values for the low-energy constants. We choose g = 0.6 36], A# = 140 MeV, 
and / = 130 MeV. The ST/ (3) splitting of the ground state D mesons is ~ 100 MeV. Assuming 
the same 577(3) splitting for the doubly charm baryons we obtain <j/\ = — (337MeV)~ 3 . We 
shall take the strange sea quark to be degenerate with its valence counterpart, m r = m s , 
and fix the value at the physical strange quark mass. This leaves three parameters that can 
be varied, the valence pion mass, m^^x = m U u, the sea pion mass, = rrijj, and the 

renormalization scale, fi. While the /i dependence in the above expressions is exactly canceled 
by /i-dependent counterterms at this order, we have not included these counterterms, nor 
do we have a reliable way to estimate them. Hence we vary fi to effectively obtain a range 
for the unknown counterterm contributions. 

In Fig. we assess the non-analytic chiral corrections to the hyperfine splitting of non- 
strange S* and S baryons in PQ^PT. To see the effect of partial quenching on this observable, 
we plot the difference between the chiral non-analytic contribution to the hyperfine splitting 
in PQ^PT and the corresponding non-analytic correction in %PT: 

5 XHF = (Mi 3 J 2) - Mi 3/2) ) - (m£ /2) - MH /2 A , (27) 

as a function of m^sea. for a few different values of m nV3 x. For each value of m nvSL \, we let 
m nsea, range from m nV3 x up to the mass of eta-strange, m Va = m ss = m rr ps 700 MeV. Thus 
the smallest value of m^-sea corresponds to an unquenched theory at = m-Kvai, while the 
largest value of m^ea corresponds to a completely degenerate sea at the strange quark mass. 
The bands correspond to varying \x from 500 MeV to 1500 MeV, which is chosen to be the same 
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for both xPT and PQ^PT. For the xPT hyperfine splitting, we use the physical Goldstone 
boson masses, so by construction 5xhf = when m,r Sea = m nvai i = m^. As demonstrated 
by Fig. chiral non-analytic corrections to the hyperfine splitting are sizable, and increase 
with increasing m nsea . Additionally corrections are enhanced as m nva [ increases. While we 
do not trust the chiral expansion at large pion masses, we have kept these in the plot to 
show that chiral corrections can become large and tend to push the hyperfine splitting down 
by tens of MeV. 

We can use our results to check the chiral corrections to the mass splitting relation 
predicted by heavy quark-diquark symmetry [H H, 0| 



M** - M= 



(28) 



In Fig. El we plot the difference between the nonanalytic chiral corrections to this relation 
in PQxPT and xPT, defined by 



$Xqdq 



M [ l /2) 



m| /2) - 1 (m p 3 / 2) 



m p 3/2) ) 



m| 3 / 2) - Mi" 



in the non-strange sector. Here ML and Mp' z> are the chiral contributions to the masses 



(3/2) 



- ( m£ /2) - M< 3/2) 



PQxPT 



X PT 



(29) 



(3/2) 



of the vector and pseudoscalar heavy mesons, respectively [32j . As in Fig. |21 we vary \i from 
500 MeV to 1500 MeV for three different values of m nvaj i. Fig.|3]shows that the chiral corrections 
to the heavy quark-diquark symmetry prediction are remarkably small compared to chiral 
corrections to the hyperfine splittings and much less sensitive to fi. It was demonstrated in 
Ref. 0| that chiral corrections to the quark-diquark symmetry prediction for the hyperfine 
splittings are small and insensitive to fi. Interestingly, this conclusion is not changed as the 
pion mass is increased, nor is it altered by partial quenching. While the hyperfine splittings 
of heavy mesons and doubly heavy baryons receive considerable chiral corrections for larger 
pion masses, the quark-diquark symmetry relation fortuitously does not. Lattice calculations 
of this quantity will not suffer large uncertainty due to chiral extrapolation. 



B. QxPT 

In the quenched theory we can similarly calculate the masses of doubly heavy baryons to 
one-loop order. The tree level contribution from 577(3) breaking counterterms is 

AfJJ = = -a Q (m q ) a . (30) 
The tree-level baryon mass splittings are 

A ba = M 3b - M Sa = M=* - M=* = a%m q ) a - (m q ) b ], (31) 

A* a = M s; - M Ea = ^A Q H + a%m q ) a - (m q ) b ). (32) 

The differences in low-energy constants notwithstanding, the one-loop mass contributions 
in QxPT are considerably altered from those in PQxPT and xPT. The sunset diagrams in 
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FIG. 3: Sxqdq as a function of m nsea for different values of m^^i- The width of the bands is the 
result of varying fi between 500 MeV and 1500 MeV. 



4* 



FIG. 4: Additional loop graphs which contribute to the masses of the doubly heavy baryons in 
QxPT. The single and double lines correspond to H baryons or H* baryons, respectively, while the 
dashed lines correspond to mesons. The filled squares denote the non-singlet axial coupling, and 
the cross denotes the axial hairpin interaction given in Eq. (jl5|) . 



Fig. ^ again contribute, however, it is clear that only hairpin contributions survive quench- 
ing. There are additional one-loop diagrams present in the quenched theory that have no 
counterparts in PQxPT or %PT. These diagrams are depicted in Fig. H]and involve the axial 
hairpin interaction in Eq. (fl3j) . Combining the hairpin diagrams and axial hairpin diagrams, 
we arrive at the one-loop quenched mass contributions 
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2g Q g 



>Q 



(47T/) 2 



5 4 
-K(0, m aa , fj) + -K(-A* aa , m aa , /i) 



3(4tt/) 2 dm 2 a 
for the spin-3/2 doubly heavy baryons, and 



5 4 
-^(0, m aa , /I) + -K(-A* aa , m aa , fj) 



MX 



(3/2) 



2g Q g' c 



(47T/) 2 



-K(0, m aa , /j,) + -K(A* aa , m aaj /i) 



3(4vr/)2 dml 
for the spin-1/2 doubly heavy baryons. 



1 8 

-K(0, m aa , fi) + -K(A* aa} m aa , fi) 



(33) 



(34) 
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FIG. 5: The quenched chiral correction to the hyperfine splitting, Mi„ — Mi , as a function 
of the quenched pion mass. The xPT band shows the corresponding range of chiral corrections to 
the hyperfine splitting in QCD. 



With the quenched chiral non-analytic corrections in hand, we can investigate their im- 
pact on the spectrum of doubly heavy baryons as we did for the partially quenched theory. 
To determine these corrections, we must estimate the parameters in the quenched chiral La- 
grangian. Little is known about these parameters. We adopt a central value of /io = 700 MeV 
from the analysis of Ref. [33|. The remaining parameters are taken to be the same values 
as in xPT, though there is no justification for this. Finally we use two values, g'® = ±0.6, 
for the parameter that has no analogue in the unquenched theory. In Fig. we plot the 
quenched chiral correction to the hyperfine splitting, — M~^ 2 \ as a function of the 

pion mass. The band plotted corresponds to a variation of \x from 500 MeV to 1500 MeV. 
As in the partially quenched theory, we see that the mass splitting receives sizable chiral 
corrections that increase with the pion mass. The trend is the same as in the partially 
quenched theory provided and g'® have the same sign: quenched lattice calculations at 
larger pion masses than that in nature will underestimate the hyperfine splitting. When 
g® and g'® have opposite signs, it is possible for the splitting to be overestimated. Lastly 
we can test how well heavy quark-diquark symmetry holds up in the presence of quenched 
chiral corrections. In Fig. El we plot the quenched chiral correction to the quark-diquark 
symmetry relation for the hyperfine splittings, 

A % DQ = Mg ,2) - M { i ,2) - 1{M$M - i\4 3/2) ) , (35) 

as a function of the pion mass, again varying fi from 500 to 1500 MeV. Surprisingly, the 
quenched chiral corrections do not sizably alter the quark-diquark symmetry prediction. 
The corrections are only a few MeV. The \i dependence of the result is minimal and in fact 
one can show analytically that the \x dependence vanishes for this observable if g'® = 0. 



IV. ELECTROMAGNETIC TRANSITIONS IN PQxPT 

As a final application we determine the electromagnetic transition moments of the doubly 
heavy baryons in PQxPT and %PT. To include electromagnetism into the theory, we merely 
gauge a background U(l) field. In doing so, recall that the light-quark electric charge matrix 
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FIG. 6: The quenched chiral correction to the heavy quark-diquark symmetry breaking parameter, 
A-Qqdq, as a function of the quenched pion mass. The band shows the corresponding range 
of chiral corrections to quark-diquark symmetry in QCD. 




FIG. 7: One-loop graphs contributing to the doubly heavy baryon electromagnetic transitions in 
PQxPT. The single lines correspond to H baryons, double lines to H* baryons. The dashed lines 
correspond to mesons, while the wiggly lines correspond to photons. The coupling of Goldstone 
mesons to baryons is the non-singlet axial coupling given in Eq. ()16|) . 



Q is not uniquely defined in PQQCD [37]. The only constraint one must impose is that the 
charge matrix Q has vanishing supertrace. Following Refs. [sH 2^1 5 we use 

Q = diag (q u , q d , q s , q j} q h q r , q u , q d , q s ) , (36) 

along with the condition qj + qi + q r = 0. QCD is recovered in the limit of degenerate valence 
and sea quarks only for the particular choice: q u — qj — |, and q d = q s = q\ = q r = — |. 
Letting the charges be arbitrary, however, enables us to track the flow of charge in loop 
diagrams. 

The relevant electromagnetic couplings in the chiral Lagrangian with quark-diquark sym- 
metry are contained in the terms [13| 

C = f (H-K Q <r ■ B) + g («•£ • B U) + ^ (H^HQS) V,,^ , (37) 

where Q' is the charge of the heavy quark and the 5x5 matrices are 

£ ~=(V-2°tJ- < 38 ) 
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The first term in Eq. ()37|) gives the leading contribution to the decay rate in the 1/tuq 
expansion. The contribution of this term at tree level is similar to the light quark magnetic 
moment contribution to the transition in the non-relativistic quark model, i.e. the quark 
model reproduces the tree-level result of HH^PT upon the identification l/m q — > (3. The 
second term in Eq. (|37|) is the contribution from the magnetic moment of the heavy quark 
which is Oil/mq). The third term in Eq. (|3T|) is also 0(l/mg) suppressed. The notation 
{ij} denotes symmetrization and tracelessness in the indices ij. In the heavy meson sector, 
this term contributes to the electric quadrupole moment of the P* but not to the P* — > P7 
transition. It does give a contribution to the decay S* — > H7, which however is suppressed 
by rather than 1/?71qA x because it does not interfere with the leading-order con- 

tribution. Explicitly the formula for the decay width is given by 



\ a 



*a7) = T 



/3Qa~ 



2m c 



+ SPa) + 



l E 2 



Qa 



+ 5(3 



E2a 



Ms. 7 ' 



(39) 



where 8(3 a and 8j3E2a represent chiral corrections to the tree-level magnetic dipole and electric 
quadrupole transition moments, respectively. The loop diagrams shown in Fig. [7| give the 
leading non-analytic chiral corrections. 

The Lagrangian in Eq. (}3Tj) is easily extended to PQxPT. Using the definition of the 
super- field 7i in Eq. (jllj) . we can expand the Lagrangian in Eq. (}3Tj) to find the doubly heavy 
baryon part and thereby determine the contribution to the transition moments at tree level. 
As this contribution only involves valence quarks, the result is the same as in [l3|. 
The coefficients j3 and 0E2 have the same value as in ;\PT which can be demonstrated by 
matching to PQxPT in the sea sector. In the heavy quark limit, the operator proportional to 
(3e2 can be neglected. When 0(l/rriQ) hyperfine splittings are included in the loop diagrams, 
there are non-analytic corrections to both j3 and /3e2- At one-loop order in PQ^PT, we must 
calculate the diagrams shown in Fig. [7| Explicitly we find 



J2 (Qb-Qa) [ dxl A 1 (x)\og 



m 



ab 



b=j 

+ 5 



- m ah R 



AxQe) 



A 2 (x) log 



m 



ab 



m ab R 



(40) 




A 2 (x)m ab 



2R[-z Lj -)h ( 41 ) 



[A 2 (x)} 2 -m ab 

In expressing the transition moments, we have made the following definitions 

A 1 (x) = -A: b + xAl a (42) 
A 2 (a;) = -A ab + xA* aa , (43) 

and used the function R(x) defined in Eq. (|23j) . The result can be deduced by replacing 
the sea charges and sea masses with the values of their valence counterparts. 
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FIG. 8: Additional loop graphs for the doubly heavy baryon electromagnetic transitions in QxPT. 
The single lines correspond to 3 baryons, double lines to 3* baryons. The dashed lines correspond 
to mesons, while the wiggly lines correspond to photons. The filled circles denote the non-singlet 
axial coupling given in Eq. Q15|) . the photon coupling is proportional to the magnetic couplings in 
Eq. (|37|). and crosses denote the quenched hairpin interaction. 

Ref. [l]|] considered the loop corrections in in the approximation where Ah = 
and cr = 0. This corresponds to taking the heavy quark limit and keeping only O(^m^) 
nonanalytic chiral corrections. In this limit the result can be inferred from the calculation 
of Ref. |40j, |4jJ for heavy mesons: 

5Pa = ~£p^ Qb ~ Q ^ mab Sp E2a = 0. (44) 

7I "-' b 

For xPT with physical parameters, using Eq. (j4Tj|) instead of Eq. (}4*4"j) has little effect on 
the transition magnetic moment, (3Q a — Q'/(2mg) + 5j3 a . The transition magnetic moment 
changes by -14% (-1%) for // = 500 MeV (1500 MeV) for a = 1 and by -1% (+3%) for /i = 
500 MeV (1500 MeV) for a = 2. We conclude that keeping only the 0(^m q ) nonanalytic chiral 
correction is an adequate approximation for estimating the electromagnetic decay width of 
the S* in QCD. Away from the physical value of the quark mass, chiral corrections become 
increasingly important and hence formulae in Eqs. ()40|) and (|41|) should be considered for 
lattice extrapolations. 

The ^-integrals can be evaluated analytically as shown in the Appendix. In actual lattice 
calculations of transition moments, the current insertion method is limited by the discrete 
lattice momenta available. The matrix elements contributing to the electromagnetic decay 
are zero for vanishing photon three-momentum q. If we wish to extract these matrix elements 
from a lattice simulation, it is necessary to know the momentum transfer dependence of 
the corresponding form factors in order to perform a model independent zero-momentum 
extrapolation. Formulae for the multipole transition form factors can be recovered from our 
expressions above by merely using the replacement 



rricf, — ► \jm\ - x(l - x)q 2 1 (45) 

where q 2 is the virtuality of the photon. 

While we do not wish to encourage quenched calculations of the transition moments, 
we shall comment briefly on the oddities of the quenched chiral calculation at one-loop 
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order. None of the diagrams in Fig. [7| survive quenching. There is, however, quark mass 
dependence at this order, but it arises from the hairpin interactions depicted in Fig. |HJ 
These diagrams are enhanced over their PQxPT and ;\PT counterparts due to the size of 
the singlet parameter /1q. This was first observed for baryon magnetic moments in Ref. 42l |. 



and applied to baryonic transition moments in Ref. [43(. Thus in the chiral limit, the Qx?T 
transition moments have chiral singularities 

o Q' 

Pmi ~ /^Mi 9 + a Mi l4 logm q + ... 

rriQ 

@E2 ~ Pm + a E2 (i 2 \ogm g + . . . , 
which should be contrasted with the chiral limit behavior in both PQxPT and %PT 

R ru rt (0) ®' -l (1/2) PZT 4- 

Pmx ~ Pmi r Pmi V m <? + • • • 

m Q 

(3e2 ~ P e\ + Pm m q log m q + . . . , 

which are non-singular. 

The inclusion of dynamical quarks thus makes a dramatic difference in the chiral behavior 
of these observables. There are two types of dynamical contributions: those from the gauge 
configurations and those from the current operator self-contraction. The latter contributions 
are notoriously difficult to determine from the lattice. The effective theory can remedy this 
situation for the electromagnetic current in three-flavor simulations. Because we require 
strQ = 0, the choice qj = qi = q r = is not excluded, i.e. one can ignore the self-contraction 
of the current. In this case, one is not determining the physical transition moments, but 
rather the physical parameters (3 and /3e2 i n the effective theory. These in turn can be used 
to make physical predictions. 



V. SUMMARY 



In this paper, we have extended the chiral Lagrangian with heavy quark-diquark symme- 
try to quenched and partially quenched theories. This extension allows for the derivation 
of extrapolation formulae necessary to connect lattice QCD data to the chiral limit. These 
formulae are easily modified to include the effects of finite volume, see Ref. ji^]. Furthermore 
light-quark discretization effects can be incorporated, but depend on the particular fermion 
discretization employed in the lattice calculation. 

As an example, we determined the quark mass dependence of doubly heavy baryons in 
both quenched and partially quenched chiral perturbation theory. We demonstrated that the 
hyperfine splitting between S* and H baryons is sensitive to chiral non-analytic corrections 
both in the quenched and partially quenched theories. Thus careful chiral extrapolation is 
required to connect lattice data to the doubly heavy baryon spectrum in the chiral regime. 
Surprisingly, we found that the heavy quark-diquark symmetry relation is rather insensitive 
to chiral corrections. Despite the possible enhancement by factors of /Zq multiplying quenched 
chiral logarithms, this conclusion is true in the quenched theory. 

We also investigated the chiral corrections to the electromagnetic transition moments. 
The width of the 5* baryons should be dominated by the decay H* — > H7. The lattice can 
be used to predict the electromagnetic decay widths and our results should be helpful for 
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chiral extrapolations in these calculations. Since the doubly heavy baryon hyperfine splitting 
is expected to be « 100 MeV, lattice pion masses will not allow on-shell intermediate states 
as one nears the chiral regime. This is in contrast to the A — > A7 transitions where the pion 
mass must be m v > 300 MeV to avoid the multi-particle continuum. Furthermore, lattice 
data can test heavy quark-diquark symmetry relations between the S* — > S7 and P* — > P7 
decays. 

Investigation of heavy mesons and doubly heavy baryons on the lattice provides a way 
to explore the implications of heavy quark-diquark symmetry. By varying the light quark 
mass, one will be able to see how well the predictions hold up against chiral corrections. 
Additionally by varying the heavy quark mass, one will be able to see in what regime this 
symmetry of the strong interaction emerges. 
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APPENDIX 



Here we give some formulae which are useful for evaluating the ^-integrals in Eqs. (j4*Uj) 
and (|41|) . A useful observation is that 



dR(x) 



x 



-R(x) - 2 . 



dx x 2 — 1 

Using this it is straightforward to verify the following anti-derivatives: 



(46) 
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x 2 — 
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x 3 




x 2 — 


1 



dx R(x) 
R{x) 



1 1 



R 2 (x) 



Ax 2 - I 
R(x) = R(x) + 2x , 

R(X) = -^—R(X)+ g 



(47) 



All integrals in Eqs. (j40|) and (J4*T|) can be evaluated using these identities. For example, for 
one of the integrals appearing in 5/3 a , we find 
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m ab 
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n ab 
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where 



G(x, y) = x 2 log(y) - — - - log 2 = 

^ o \x+v^— 1 + «e / 

x\/a; 2 — 1 j / a; — \Jx 2 — 1 + ie\ 
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